% Prints out some steady-state statistics
fprintf('\n')
fprintf('Model                            : %0.1g \n',model)  
fprintf('\n')
fprintf('Annual inflation                 : %0.8g%% \n',100*(mu^mfreq-1))  
fprintf('\n')
fprintf('Frequency of price changes       : %0.8g%%  \n',[freqpchanges])
fprintf('\n')
fprintf('Fraction of price increases      : %0.8g  \n',[fracPriceIncr])
%fprintf('\n')
%fprintf('Fraction of price changes <=5%%   : %0.3g   %0.3g\n',[fracSmallChanges PctLessthan5])
%fprintf('Fraction of price changes <=2.5%% : %0.3g   %0.3g \n',[fracSmallChanges25 PctLessthan25])
fprintf('\n')
fprintf('Mean price change                : %0.8g   \n',EPchange)     
fprintf('Std of price changes             : %0.8g   \n',StdPchanges)  
fprintf('Mean abs of price changes        : %0.8g   \n',MeanAbsPchange)  
fprintf('Kurtosis of price changes        : %0.8g   \n',KurtPchanges)                                           
fprintf('\n')
fprintf('\n')
fprintf('Menu cost (kappa)                      : %0.8g   \n',par_.menucost)     
fprintf('Std of idiosyncratic shocks (std_eA)   : %0.8g   \n',par_.std_eA)  
fprintf('Information cost (alpha)               : %0.8g   \n',par_.alpha)  
fprintf('Hazard slope (alphapos)                : %0.8g   \n',par_.alphapos)  
fprintf('Labor disutility coefficient (chi)     : %0.8g   \n',par_.chi)
fprintf('Real wages (wbar)                      : %0.8g   \n',par_.wbar)
fprintf('Hazard minimum (omega)                 : %0.8g   \n',par_.omega)
fprintf('\n')
fprintf('\n')
fprintf("Intensive (impact,relative)                   : %3.2f%% (%3.2f%%)\n",intensive(1,1)*100,relIntensive(1,1)*100);
fprintf("Gross extensive (impact,relative)             : %3.2f%% (%3.2f%%) \n",grossExtensive(1,1)*100,relGrossExtensive(1,1)*100);
fprintf("Selection (impact,relative)                   : %3.2f%% (%3.2f%%) \n",selection(1,1)*100,relSelection(1,1)*100);

fprintf("Intensive (1m,relative)                   : %3.2f%% (%3.2f%%)\n",intensive(1,2)*100,relIntensive(1,2)*100);
fprintf("Gross extensive (1m,relative)             : %3.2f%% (%3.2f%%) \n",grossExtensive(1,2)*100,relGrossExtensive(1,2)*100);
fprintf("Selection (1m,relative)                   : %3.2f%% (%3.2f%%) \n",selection(1,2)*100,relSelection(1,2)*100);


fprintf('\n')
fprintf("Intensive (12m,relative)                   : %3.2f%% (%3.2f%%)\n",intensive(1,12)*100,relIntensive(1,12)*100);
fprintf("Gross extensive (12m,relative)             : %3.2f%% (%3.2f%%) \n",grossExtensive(1,12)*100,relGrossExtensive(1,12)*100);
fprintf("Selection (12m,relative)                   : %3.2f%% (%3.2f%%) \n",selection(1,12)*100,relSelection(1,12)*100);

fprintf('\n')
fprintf("Intensive (24m,relative)                   : %3.2f%% (%3.2f%%)\n",intensive(1,24)*100,relIntensive(1,24)*100);
fprintf("Gross extensive (24m,relative)             : %3.2f%% (%3.2f%%) \n",grossExtensive(1,24)*100,relGrossExtensive(1,24)*100);
fprintf("Selection (24m,relative)                   : %3.2f%% (%3.2f%%) \n",selection(1,24)*100,relSelection(1,24)*100);

if model==3                                                        
   fprintf('\n')
   fprintf('Menu cost as %% of revenue        : %0.3g%% \n',McostinRev)
   fprintf('Menu cost as %% of wage bill      : %0.3g%% \n',McostinWageBill)
   fprintf('Menu cost as %% of consumption    : %0.3g%% \n',McostinCbar)
end
fprintf('\n')